Collective many-body interaction in Rydberg dressed atoms 
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We present a method to control the shape and character of the interaction potential between 
cold atomic gases by weakly dressing the atomic ground state with a Rydberg level. For increasing 
particle densities, a crossover takes place from a two-particle interaction into a collective many-body 
interaction, where the dipole-dipole/van der Waals Blockade phenomenon between the Rydberg 
levels plays a dominant role. We study the influence of these collective interaction potential on a 
Bose-Einstein condensate, and present the optimal parameters for its experimental detection. 



The remarkable success of ultra-cold atomic gases in 
exploring quantum phenomena from the weakly inter- 
acting regime to strongly correlated many-body physics 
is founded in the microscopic understanding of the two- 
body interaction potential, and the possibility to control 
it by external fields. Several tools are nowadays avail- 
able to control the strength and shape of the interaction 
potential with the most prominent examples beeing mag- 
netic and optical Feshbach resonancesllj-y] , atomic gases 
with large magnetic dipole moments [5| , as well as the re- 
cently realized polar molecules [6] . An alternative class of 
strong interactions appears between atoms excited into 
a Rydberg state Q- In this letter, we show that the 
shape and character of the interaction potential between 
ground state atoms can be altered by dressing them with 
a Rydberg level giving rise to a collective many-body in- 
teraction. 

The optical dressing of ground state atoms with an 
excited electronic state has extensively been studied in 
the past in the context of "blue shield" for the supres- 
sion of inelastic collisions [§]. Motivated by recent ex- 
perimental progress in manipulating Rydberg atoms also 
in the ultra-cold regime j9i] , the possibility of controlling 
the interaction via Rydberg dressed ground state atoms 
has been proposed for the dipole-dipole interaction in the 
strongly interacting regime [13], and for inducing short 
range interactions in a Bose-Einstein condensate . In 
contrast to the earlier proposals [12[, these approaches 
use a finite detuning in order to reduce losses via spon- 
taneous emission from the excited Rydberg state. 

Here, we study the influence of Rydberg dressed in- 
teractions on the ground state wave function of a Bose- 
Einstein condensate. The combination of strong interac- 
tions between the Rydberg state and the detuning gives 
rise to a large effective range of the interaction poten- 
tial, which exceeds the natural interparticle distance in 
cold atomic gases. Within this regime, we show that a 
crossover from a two-particle interaction potential to a 
collective many-body interaction takes place. This col- 
lective many-body interaction appears as a consequence 
of the dipole/van der Waals blockade extensively studied 
in the past [l3l-[l7|. We provide a description for this 
collective interaction in a Bose-Einstein condensate and 



discuss its experimental signatures. 
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FIG. 1. (a) Setup for the ground state dressing with the Ry- 
dberg level \r) via the intermediate \p) state with the cor- 
responding Rabi frequencies and detunings. (b) Diagram 
showing the crossover from the two-body interaction poten- 
tial into the regime with a collective many-body interaction. 
The crossover line is determined by A ~ — 0*^^ (red dashed 
line), which preempts the condition of weak dressing Q, < |A|. 

In the following, we are interested in a two-photon cou- 
pling of the ground state atoms to a high lying Rydberg 
state; see for such an experimental setup. Then, the 
relevant internal structure for each atom is given by a 
three- level system: the atomic ground state \g)i is cou- 
pled to the Rydberg level \r)i via the intermediate \p)i 
state, see Fig. [TJ The Rabi frequency and detuning for 
the transition from the ground state to the p-level is de- 
noted as rip and Ap, and for the transition from the 
p-level to the Rydberg state fi^ and A^, respectively. 
The intermediate p-level is far detuned, Ap » f2r,rip, 
and can be adiabatically eliminated. Then the dynam- 
ics for a single atom reduces to an effective two level 



system Hi — h ^a, 
Rabi frequency 
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/2 with the two-photon 



f7ri^p/2|Ap| and the total detuning 

A = Ar -I- Ap; here tXa'' denote the Pauli matrices. In 
the following, we first derive the two-body interaction po- 
tential in the regime with large detuning, i.e., $7 <^ |A|. 
Then, the ground state atom is weakly dressed with the 
excited Rydberg level \+) « (1 - Vt^ /d,J\^)\g) + 0/2A|r) 
providing a small fraction / w ri^/4A^ of excited Ry- 
dberg atoms. At large distances and in absence of 
Forster resonances 15], the interaction between Ryd- 
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berg states is dominated by the van der Waals interac- 
tion K-dw(r) = C^/r^. In the following, we are inter- 
ested in a repulsive van der Waals interaction Cg > 
and red detuning for the two-photon transition A < 0; 
the generalization to attractive interactions with blue de- 
tuning is straightforward. Then, the weak dressing of 
the two-ground state atoms gives rise to the induced 
interaction V^tt(r) ^ fi^'/ieA^'C'e/'r^ at large distances 
r > ^0 = (C6/2ft|A|)^/^ On the other hand, for short 
distances r < the two atoms enter the van-der Waals 
blockaded regime, where only one atom is dressed with 
the Rydberg state, and the interaction saturates at an en- 
ergy shift ^ fi^'/A^Ce/^Q. The full interaction potential 
Kitf is derived within the adiabatic Born Oppenheimer 
approach K„(r) = hO." / + [r/^of]'^ [IS El, 
which in Fourier space reduces to 



K.«(q) 
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with f{z) = [e-^/2-2cos(V3z/2 + 7r/3)] /z. 

Within first Born approximation, this two-body inter- 
action gives rise to the s-wave scattering length 
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Note, that on short distances, additional crossings with 
the surrounding Rydberg levels appear, which give rise to 
very narrow resonances and opens up an additional loss 
chanel. However, due to the extremely narrow struc- 
ture of these resonances one can expect that the par- 
ticles mainly move diabatically across, and these losses 
are strongly suppressed; a detailed theoretical analysis of 
these losses is so-far missing. 

In the many-body system, the effective interaction be- 
tween the particles can in general also contain many- 
body interaction terms. For weak Rydberg dressing with 
(ri/A)^ ^ 1, these additional terms become relevant at 
high particle densities n: first, the two-particle block- 
ade radius gives rise to a maximal density of excited 
Rydberg atoms l/^g above which collective phenomena 
become relevant. On the other hand, the density of ex- 
cited Rydberg atoms due to the weak dressing is given by 
by (ri/2A)^n. Consequently, the validity of the two-body 
interaction Kff(r) is limited to the dilute regime 
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At higher densities n > = / {Vi^ ^/C^) 

a crossover into the collective regime takes place (see 
dashed line in Fig.[T]) and the effective interaction is dom- 
inated by collective many-body interactions. The behav- 
ior I A|^ ^ (Cen^ri'') of this cross-over line derived within 
this simple estimate agrees with the one expected from 
the universal scaling theory 3 19| . 



In the following, we provide a description of this ef- 
fective interaction for a Bose-Einstein condensate. The 
idea is to derive for a fixed atomic density n the inter- 
nal energy E^^i [n] of the driven Rydberg system within a 
Born- Oppenheimer approximation. This internal energy 
then describes the collective interaction potential giving 
rise to a generalized Gross-Pitaevskii equation 

ihdt^{t, r) = {Ho + gsn{t, r) + [n{t, r)]} r) (4) 

with the condensate wave function ip{t,r), the density 
n{t,r) = |^(t,r)p, and the derivative i?'fjn] — dnE^ff[n]. 
Here, Hq = — ?i^A/(2m) -I- K,ct(r) — ^ describes the non- 
interacting part with V^xt(r) the external trapping po- 
tential, m the mass of the particles, and fi the chemical 
potential, while gs accounts for the s-wave scattering be- 
tween ground state atoms. This approach remains valid 
for momenta of the wave function smaller than the char- 
acteristic length scale for the effective interaction. 

The internal energy E^f{[n\ can be derived within a 
mean-field/variational approach. The Hamiltonian for 
the internal degree of freedom is described by a spin 
Hamiltonian uM 



H 



2 tz" l"^' 



with p(*^ = (1 -I- ai^^)/2 the projection onto the ex- 
cited Rydberg state \r). The strong repulsion between 
the Rydberg atoms gives rise to the correlation function 
g{Yi — Yj) = (P(')p*^-'))//^. The correlation function van- 
ishes on short distances, i.e., g{r) « for |r| <^ ^, due 
to the Blockade of Rydberg excitations, and approaches 
unity g{r) — >• 1 at large distances |r| ^. The charac- 
teristic correlation length scale ^ denotes the collective 
blockade radius. Note, that this correlation function is 
independent on the remaining particle positions due to 
the homogeneous distribution of the particles in the BEG. 
A variational wave function, which allows one to include 
the strong correlations takes the form 
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where Css'(r) with s e {ff, ^} account for correlations 
on distances shorter than a characteristic length scale 
^, and reduce to Css' — > ctgOis' for large distances with 
ag — cos 9 and = sinO. Then, this wave function 
describes a paramagnetic phase with the spins aligned 



and 



along the direction (cr^*-'), i.e., cos26' = — (ct, 
sin26' = — (fTx''), and the fraction of excited Rydberg 
atoms reduces to / = (P^*^) = sin^ 9. Note, that the re- 
lation between Css' (r) and the correlation function g(r) 
is highly non-trivial. In the following, we are interested 
in using the correlation length ^ as a variational param- 
eter, and therfore, the precise shape of Css' is irrelevant. 



3 



The energy density e^^^{9,^) 
degree of freedom reduces to 



|_ff|0) for the internal 



— — sm2t^ + Asin^ . (6) 



The variational parameters are the excited Rydberg frac- 
tion / = sin^ 9 and the correlation length 4. The di- 
mensionless parameter X/^^ = J d'^xg{x) — l/^"^ • 
/ cPxg{£,x.)/\x.\^ is determined by the details of the cor- 
relation function. It is important to note, that the qual- 
itative behavior of the energy functional is independent 
on the precise value of A, and therefore, this parameter 
will be fixed by the constraint to reproduce the correct 
asymptotical behavior at low densities. The dressing of 
the ground state atoms with the excited Rydberg level 
leads to an energy shift eg = ~hny^^ + TP /2 even in 
absence of Rydberg interactions. Subtracting this con- 
tribution, the energy functional £'„ff[n], describing the 
Rydberg dressed interactions in the collective regime, is 
obtained by the energy for the state adiabatically con- 
nected to the non-interacting system: this wave function 
is obtained by minimizing the energy function e„„, i.e.. 



E„ft \ri\ = min 

m} 
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for fixed ground state density n, Rabi frequency and 
detuning A. The value 6 for the minimal energy density 
follows from the equation dee^^^{0,£,) — 0; note, that this 
condition is equivalent to the self-consistency equation 
derived in Ref. On the other hand, the last term 

in equation ^ always gives a positive contribution to 
the energy functional, and it will be minimized by maxi- 
mizing the correlation length 4. However, the correlation 
length satisfies several constaints: in the weakly inter- 
acting regime the correlations due to the Rydberg exci- 
tations are limited by the two-particle Blockade radius 
[1^ . Entering the collective regime the distance between 
the Rydberg atoms is further reduced and the correla- 
tions are bounded by {l/nfY^^, i.e. S, « (l/rt/)^/"^, see 
inset in Fig. [5J Note, that the qualitative behavior of the 
system is independent on the precise choice of ^ but only 
weakly renormalized the dimensionless parameter A. 

The energy functional -B^tfl'^] is shown in Fig. 2 for 
varying density n. In the low density regime, n <C ric = 
4|A|^/^/(17^-^C6/2), it reduces to a pure two-particle in- 
teraction E^tt = g^ttTi^ /2. The comparison with the exact 
two-body interaction Eq. ([2]) allows us to fix the dimen- 
sionless parameter A and within first Born approximation 
it takes the form A = 27r^/3. Strong deviation from the 
two-body interaction can already be seen at densities as 
low as 0.05nc, see Fig. [2j Collective blockade phenom- 
ena give rise to a very broad crossover, with the energy 
functional i?ett [n] eventually saturating at a chemical po- 
tential E,ii[n] « = nh (VA^ + Q"^ - | A|) /2 at high 
densities n ^ ric- This behavior can be understood from 
the following argument: (i) each atom within a Blockade 
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FIG. 2. (a) Energy funcitonal d„E,ii[n] for |!^/A| = 0.072. 
The dashed line describes the behavior at low densities 
i5off[n] ~ (?off?i^/2 and its deviation indicates the early break 
down of the two-particle interactions, while the dotted line 
accounts for the saturation at high densities, although this 
saturation is very slow. The inset shows the behavior in the 
Blockade radius across the critical density ric- (dotted line 
is guide to the eye), (b) Behavior of the density profile for 
a Bose-Einstein condensate with Qs = and the collective 
many-body interaction E'^n [n] for increasing particle numbers 
in a harmonic trap with lj = 20Hz. A characteristic feature 
is the strong increase of the central density n(0) in contrast 
to the conventional Thomas- Fermi scaling n(0) ~ N^''^. 



radius gives rise to an energy shift in the functional, and 
(ii) in the high density limit to fraction of excited Ryd- 
berg atoms vanishes. Therefore, all atoms up to a non- 
extensive part are within a Blockade radius, and thus 
basically free. 

Finally, we can derive the ground state wave func- 
tion for a Bose-Einstein condensate in an harmonic trap 
with the infiuence of the Rydberg dressed interaction 
within the Thomas-Fermi approximation. The density 
profile for the collective many-body interaction is shown 
in Fig. Hfb) for increasing particle number A'^. A charac- 
teristic feature is the strong increase of the density in the 
trap center. This behavior is a consequence of the linear 
behavior of the collective interaction potential at high 
particle densities n 3> Uc, and allows one to distinguish 
the collective interaction from a two-body interaction or 
heating effects. The density profile for an atomic conden- 
sate of ®^Rb atoms with a background scattering length 
as = 5.7nm and fixed particle number N — 10^ is shown 
in Fig|3l For weak Rydberg dressed interactions, the 
density profile is given by the inverted parabola profile 
n{r) = — Kxt(i")] /ffcff- The crossover into the collec- 
tive interaction regime is again signaled by the charac- 
teristic change in the shape of the atomic cloud. 

In the following, we will determine the parameters for 
an experimentally realistic setup. For weak dressing, the 
population of the Rydberg state is suppressed by the fac- 
tor (ri/A)^, as well as the intermediate p-level by the 
factor (Op/Ap)^. Denoting the decay rate of the Ryd- 
berg state and the p-level with 7^ and 7^, respectively, 
the total decay rate 7„ff for spontaneous emission is given 
by 
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FIG. 3. Density profile n(r) for a Bose-Einstein condesate of 
*''Rb atoms in a fiarmonic trap uj = lOOHz witfi background 
scattering lengtli as — 5.7nm for increasing influence of the 
Rydberg dressed interaction. The inset shows the path within 
the parameter space for increasing interaction strengths. 



The detuning Ap and Rabi frequency flp of the laser cou- 
pling the ground state and the intermediate p-level only 
determine the global life time via its ratio (ftr/Ar)'^. 
Second, the Rabi frequency fir for the coupling to the 
Rydberg state is limited by experimental restrictions on 
available laser power. Consequently, the strongest influ- 
ence of the Rydberg dressed interaction is obtained by 
the choice of the detuning A, which extremizes the in- 
duced s-wave scattering length for fixed life time of the 
system, i.e., d^g^a = 0. Then, the optimal detuning sat- 
isfies Ac = — |r2r|\/7r/(287p) with the induced s-wave 
scattering length 
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For *^Rb weakly dressed with the Rydberg level |35s), 
the van der Waals interaction takes the strength Cg = 
1.31 lO^^a.u. [l3| , and the decay rates 7^ = 4kHz and 
7p = 6MHz. For a realistic experimental setup, Rabi fre- 
quencies with f2r — 22MHz have been reached. Finally, 
a life time of 7 = 6Hz is sufficient for a fast experiment, 
and we obtain 



1.36 7^ nmHz ^ = 49.5nm, 



(9) 



compared to = 5.7nm, which is an increase of the 
scattering length by an order of magnitude. The optimal 
detuning reduces to |Ac| = 107kHz with the character- 
istic length scale ^0 = 3/xm. The coupling to the p-level 
requires the strength |fip/Ap| — — 7.0 10^^, 

and the effective Rabi frequency becomes £7 = 7.8kHz, 
thus |f7/Ac| = 1/77/27^ = 0.072. The critical den- 
sity for the crossover from the two-body interaction to 
the collective many-body interaction reduces to ric — 
2.57 lO-'^'^cm""^, and consequently, for typical peak densi- 
ties ripeak ~ lO^^'cm"'^ of a typical BEC experiment the 



system is well within the collective regime, and allows 
for the detection of this novel interaction. Note, that 
within the collective regime, the number of excited Ryd- 
berg atoms is further reduced due to the Blockade phe- 
nomena. This in principle allows one to further reduce 
the detuning without increasing the losses from sponta- 
neous emission. However, in a trapped system, this will 
increase the losses at the edge of the condensate, where 
the density drops below the critical density ric- 



Finally, we would like to point out, that the weakly 
dressed Rydberg interaction could be further tune d by 
a micro- wave field coupling different Rydberg levels [20[ ■ 
As opposed to the much stronger van der Waals repul- 
sion between the Rydberg levels, the micro-wave field 
only gives rise to a weak dipole-dipole interaction, which 
becomes relevant on large distances r ^ ^ due to its 
slow decay. As a consequence, it does not change the 
behavior on distances comparable to the Blockade ra- 
dius. Such a scenario is in analogy to the possibilit y t o 
tune interaction potentials with cold polar molecules t2l|, 
and will allow one to explore the interesting properties 
of anisotropic dipole-dipole interactions in the collective 
regime. 



Support from the DFG (Deutsche Forschungsgemein- 
schaft) within SFB/TRR21 is acknowledged. 



[1] 
[2: 

[3: 
[6: 
[s: 

[9 

[10: 
[11 

[12 
[13 
[14 
[15 
[16 
[17 
[18 

[19 
[20 
[21 



E. Tiesinga, B. J. Verhaar, and H. T. C. Stoof, Phys. 
Rev. A 47, 4114 (1993). 

M. R. Andrews et al., Nature 392, 151 (1998). 

R Courteille et al., Phys. Rev. Lett. 81, 69 (1998). 

P. O. Fedichev et al., Phys. Rev. Lett. 77, 2913 (1996). 

T. Lahaye et al.. Nature 448, 672 (2007). 

K.-K. Ni et al.. Science 322, 231 (2008). 

K. Singer, J. Stanojevic, M. Weidemiiller, and R. Cote, 

J. Phys. B 38, S295 (2005). 

J. Weiner et al., Rev. Mod. Phys. 71, 1 (1999). 

R. Heidemann et al., Phys. Rev. Lett. 100, 033601 

(2008). 

G. Pupillo et al., Phys. Rev. Lett. 104, 223002 (2010). 
N. Henkel, R. Nath, and T. Pohl, Phys. Rev. Lett. 104, 
195302 (2010). 

L. Santos et al., Phys. Rev. Lett. 85, 1791 (2000). 
K. Singer et al., Phys. Rev. Lett. 93, 163001 (2004). 
D. Tong et al., Phys. Rev. Lett. 93, 063001 (2004). 
T. Vogt et al., Phys. Rev. Lett. 97, 083003 (2006). 
R. Heidemann et al., Phys. Rev. Lett. 99, 163601 (2007). 
U. Raitzsch et al., Phys. Rev. Lett. 100, 013002 (2008). 

H. Weimer, R. Low, T. Pfau, and H. P. Biichler, Phys. 
Rev. Lett. 101, 250601 (2008). 

R. Low et al., Phys. Rev. A 80, 033422 (2009). 

M. Miiller et al., New J. Phys. 10, 093009 (2008) 

A. V. Gorshkov et al., Phys. Rev. Lett. 101, 073201 

(2008) 



